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Opinion dynamics of random-walking agents on finite two-dimensional lattices is studied. In the 
model, the opinion is continuous, and both the lattice and the opinion can be either periodic or 
non-periodic. At each time step, all agents move randomly on the lattice, and update their opinions 
based on those of neighbors with whom the differences of opinions are not greater than a given 
threshold. Due to the effect of repeated averaging, opinions first converge locally, and eventually 
reach steady states. Like other models with bounded confidence, steady states in general are those 
with one or more opinion groups, in which all agents have the same opinion. When both the lattice 
and the opinion are periodic, however, metastable states, in which the whole spectrum of location- 
dependent opinions can coexist, can emerge. This result shows that, when a set of continuous 
opinions forms a structure like a circle, other that a typically-used linear opinions, rich dynamic 
behavior can arise. When there are geographical restrictions in reality, a complete consensus is rarely 
reached, and metastable states here can be one of the explanations for these situations, especially 
when opinions are not linear. 

PACS numbers: 89.75Fb, 87.23. Ge, 02.50.Ey, 05.40.Fb 



I. INTRODUCTION 

The attempt to investigate social systems by physicists 
is several decades old [J, even though social dynamics 
has become a popular subject in statistical physics only 
recently [2]. One of the reasons for this interest is that 
one basic approach to study social systems is similar to 
what statistical physicists typically try to do: namely, 
finding macroscopic behavior or emergence from dynam- 
ics of microscopic entities [3J 0]. While physical systems 
deal with particles, entities that make up social systems 
are humans, or groups of humans. Figuring out dynamic 
behavior of even one human being is not an easy task, but 
some aspects of collective behavior of many individuals 
are known to be describable using microscopic models [3J, 
and even universal[5]. In addition, due to the current 
ubiquity of the Internet, especially the popularity of so- 
cial networks, and the increased capability of processing 
vast amount of social data, this kind of approach has 
become not only possible, but also useful. 

Opinion dynamics is one of the social-dynamics prob- 
lems that can be closely related to physical problems. 
Microscopic models we are interested in here typically 
evolve with discrete time steps, and have the fixed num- 
ber of "agents" (actors, or individuals) with their own 
opinions. We can categorize these models using sev- 
eral basic features. Opinions can be discrete [6HT6], or 
continuous [T21 [T7rI2"5] . Examples of discrete opinions are 
yes or no on a question (2 values), evaluation on a scale 
from 1 to 5 (5 values), choices in elections (2 or more 
values), and so on. When there are more than a few 
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choices, however, continuous opinions can be used: fine- 
scaled evaluation on something on a scale from to 1, 
political views, and so on. An opinion can be a vec- 
tor of integers[THl I21)l - f2"8] as well. Another important 
feature is how a model restricts interacting partners of 
an agent at a given time. Agents typically have ongo- 
ing relationships with others, and interact with selected 
peers out of related ones. The structure of these re- 
lationships plays an important role in social dynamics, 
and networks can be used to describe these relations (so- 
cial networks). We can divide models into three dif- 
ferent cases: (i) fully-connected networks, where each 
agent is related to all other agent at any moment (there 
is no restriction, and the network concept is not neces- 
sary) [I21ITgH2rJll27H2g]: (ii) fixed networks, where each 
agent is related to the limited number of agents given 
by time-independent networks [§£□] 01H2U [Ml HZ]; 
(iii) evolving networks, where network structures evolve 
with time HM3 HI HI 123 W\- Finally, models can 
be differentiated by how updating agents are chosen at 
each time step. One can update one agent at a time 
(the serial update), or all agents synchronously, espe- 
cially when the order of update doesn't play a role (the 
parallel update) [TH1 [3T] . 

The model introduced here uses continuous opinions 
with evolving networks, and the parallel update. Agents 
reside and move randomly on a two-dimensional (2D) 
lattice. At each time step, agents update their loca- 
tions in the lattice using the 2D random walk, and 
change their opinions synchronously. Only nearest- 
neighbor interactions are allowed for opinion changes; 
hence interacting partners can be represented by a con- 
tact network\31\. whose evolution is only governed by 
the movements of agents. This is basically the process 
of repeated averaging[T71 and opinions have a ten- 
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dency to move toward those of neighbors. Our model also 
uses the threshold to restrict interactions between agents 
with the big difference of opinions (bounded confidence, 
and the use can be justified in many aspects[331 I3~4]). 
In most models with the above setup [TTJ [T^l IT81 - 
EU EH EH EH ES] , the system eventually reaches one of 
steady states, where one or more groups of agents reach 
their consensus. Unlike other models, we assume that 
both the lattice and the opinion can be periodic. The 
shape of the lattice can be either rectangular or toroidal: 
two of the simplest shapes, and yet different topologi- 
cally. Opinions can be periodic, too, when an opinion is 
about a periodic subject like the time of the year. When 
both the lattice and the opinions are periodic, we observe 
some periodic metastable states. In these states, there is 
no consensus even though opinions converge locally, and 
the whole spectrum of opinions, which depend only on 
spatial locations of agents, can coexist. The main pur- 
pose of this work is to 

The plan of this paper is as follows. In Sec. II, the 
model is described, and steady states of an extreme case 
of an one-dimensional lattice are found. In Sec. Ill, nu- 
merical results when both the lattice and the opinion are 
periodic are shown. Finally, in Sec. IV, possible exten- 
sions, and many aspects of this model are discussed. 

II. MODEL 

We propose a simple model for opinion dynamics of 
random-walking agents with only nearest-neighbor inter- 
actions. We assume that there are N agents, each with 
an opinion, and that they reside on finite 2D lattices. 
Opinion changes will come only from interactions with 
neighbors. Time is discrete, and is represented by a di- 
mensionless quantity, t, which is a non-negative integer. 
Agent i (1 < i < N) will have a location and an opinion 
at time, t. Because the structure of the lattice can play 
an important role in the dynamics, we consider two struc- 
tures: a rectangle (non-periodic) and a 2D torus (peri- 
odic) [see Fig.HTa)]. The location on the lattice for agent 
i at time t is (Xi,yi), where X{ and yi are non- negative 
integers (0 < x, < X and < < Y). 

In our model, an opinion of agent i, (pi, is a real number 
between and 1 (<fo € [0, 1]), and can be either periodic 
or non-periodic [see Fig.[ljb)]. For non-periodic opinions, 
<f>i is a number on a line between and 1; while, for 
periodic opinions, 4>i is a point on a circle (for this case, 

and 1 are the same opinion) . Then, the state of agent 

1 at time t, Si(t), is represented by three numbers, 

s i (t) = {x i (t),y i (t); ( j> i (t)}, (1) 
and the state of N agents is an iV-tuple of Sj's, 

s = (si,s 2 , ...,s N ). (2) 




FIG. 1. Structures of the lattice and the opinion. Note that, 
for example, the differences of opinions between 0.2 and 0.8 
are different for two cases: 0.6 and 0.4, respectively. 



are two kinds of movements: changing locations in the 
lattice, and change of opinions. First, we look at how 
agents move in lattices. For simplicity, we use an inde- 
pendent random- walk motion for each agent. The move- 
ment of agent i (1 < i < N) is represented by possible 
choices of locations at a next time step as follows, 

{{Xi + l,yi) 
(x i ,y l + 1) 
(ii - 1, . (3) 
(xhVi - 1) 

Here an agent can move to five neighboring locations, in- 
cluding an option for staying[36 with equal probabilities, 
1/5. If an agent has less than five choices to move, it will 
move to one of those locations with equal probabilities as 
well. For the rectangular lattice, for example, if an agent 
is on an edge, it has four possible choices with equal 
probabilities, 1/4. We assume that more than one agent 
can reside in the same location, and we define neighbors 
as agents residing at the same location at a given time. 
Since agents move randomly, neighbors of an agent will 
also change with time. This contact network consists of 
disconnected cliques of various sizes as in Ref . 1311 

Interacting partners of an agent are further reduced by 
the use of the threshold d as in other models of bounded 
confidence [13 E3 [EMU El El ESI 133 . If the difference 
of opinions for a given pair of agents is greater than d 
there will be no influence. In other words, if we introduce 
A.ji as the difference of opinions between agents i and j, 
agent i will not be influenced by agent j if \Aji\ is greater 
than d. For non-periodic opinions, Aji can be obtained 
by subtraction, 



How do states of agents evolve? The change of Si will 
depend on states of other agents. For each agent, there 



Aji = cpj - (j),. (4) 
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For periodic opinions, on the other hand, we set 



(a) 



t=0 



1 (if fa -fa < -0.5), 
1 {if fa -fa >0.5), (5) 
(otherwise), 



[see Fig. [TJb)] . The ranges of both A 3 » and (i are differ- 
ent for two types of opinions: for non-periodic opinions, 
— 1 < Aji < 1 and < d < 1, and for periodic opinions, 
-0.5 < Aj, < 0.5 and < d < 0.5 (d = 0.5 for peri- 
odic opinions is excluded because an uncertainty can be 
introduced.) 

Then, we can define the set of neighbors with similar 
opinions as Ki, where agents in this set are neighbors of 
agent i at a given time and absolute values of differences 
of their opinions with agent i are less than or equal to d. 
By including agent i in K{ (i G JQ), there will be at least 
one element in i-Q. Then, the opinion of agent i, fa, at 
t + 1 becomes 



fa(t + l) =fa{t)+g 



E 



jeKi(t) AjiM 



(6) 



where g is a convergence parameter (0 < g < 1), and 
|-Ki(t)| is the number of elements in the set Ki at time 
t. When \Ki\ = 1, there is no interaction for agent i, 
and fa will not change. When \Ki\ = 2, the interaction 
is binary, and Eq. |6) becomes the equation in the Def- 
fuant model 29J (in the original model, y, has been used, 
which is basically g/2). When \Ki\ > 2, agent i is inter- 
acting with more than one neighbors at once as in the 
Hegselmann-Krause model[T8]. 

The right-hand side of Eq. (|6| can be also written as 
fa + g(fa — fa), where fa is the average opinion of agents 
in Ki. When g — 1 and d is maximal, all interacting 
agents will have the same averaged opinion at the next 
step. For example, for non-periodic opinions, when two 
agents with fa and fa interact, both of their opinions will 
become (fa +fa)/2 at the next step. Note that care has 
to be taken when the opinion is periodic. For example, 
the average of opinions, 0.1 and 0.9, should be 0, not 0.5 
as in the non-periodic case. For periodic opinions, the 
value of fa can become greater than 1 or less than after 
Eq. ([6| is applied; in those cases, we can adjust fa value 
by subtracting or adding 1 to keep fa in [0,1]. When 
< g < 1, every interacting agent will move toward a 
certain value at the next step. If the same agents interact 
for more than one time step, their opinions will gradually 
converge to one opinion value. The bigger the value of 
g, the faster they will converge. We are not considering 
the case with g = or d = 0, because the model becomes 
trivial. This model has five parameters: N, X, Y , d, and 

a- 

At t = 0, the state of all agents is given, and opinions 
and locations thereafter will be calculated from the state 
of the previous time step. At each time step, agents un- 
dergo random walks on the lattice, and change opinions 
according to Eq. ^ synchronously. This dynamic pro- 
cess is stochastic because of random walks, even though 
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FIG. 2. Scatter plots to show time evolutions in ID and 2D 
cases with the non-periodic lattice and periodic opinions, (a) 
ID case: N = 1000, X = 1000, Y = 1, g = 1. d= 0.49; (b) 
2D case: N = 10000, X = 100, Y = 100, g = 1, and d = 0.49. 



opinion changes are deterministic. This also is a Markov 
process because the state at the previous time step is all 
we need to find the current state and beyond. 

Before looking at numerical results, we can get some in- 
sights by looking at an extreme case: the one-dimensional 
(ID) case (by setting Y as 1) with g = 1 and the max- 
imal d. In this case, opinions of agents at one location 
converge to the same value at each time step according 
to Eq. and due to local interactions, the whole state 
can be approximated by a continuous ID curve in (x, fa- 
space [see Fig.^a)]. Then, we can focus only on dynam- 
ics of this curve instead of s in Eq. (2) . The local stability 
at a location can be achieved if the curve is locally linear 
(see Appendix for detail), and we get the global stability 
when the curve is linear everywhere, which means that 
steady states are straight lines in (x, fa). Given an initial 
condition, the system will eventually reach one of states 
represented by straight lines (it can be seen as a straight- 
ening process of a curve), and opinion values with respect 
to x will not change. The most common is the one where 
every agent has the same opinion (complete consensus), 
as most models of opinion dynamics have found. Once 
reached, opinions of all agents will not change afterward. 
We will call these steady states flat here, because they 
can be represented by flat lines in (x,fa. If we use the 
term an "opinion group" (a cluster, or a party) as a set 
of agents that have reached a consensus, there will be 
only one opinion group for maximal d, while there can 
be more than one opinion groups when d is small (say, 
d < 0.3), as seen in models of bounded confidence. 

When both the lattice and the opinion are periodic, 
however, non-flat steady states can appear. They are 
called non-flat, because they can be represented in (x,fa 
as lines with non-zero slopes. Unlike flat steady states, 
lifetimes of these states are finite in general depending 
on some parameters (will be discussed in detail later); 
therefore, these states will be also called metastable states 
(see Ref . [TBI for another type of metastable states). Note 
that non-flat states cannot be sustained in non-periodic 
cases due to the boundary effect [37;. 

It is not hard to generalize the ID results to 2D lattices, 
and steady states will be represented by 2D planes in 



FIG. 3. Scatter plots for the case with ID periodic lattice 
and periodic opinions (TV = 1000, X = 1000, Y = 1, g = 1, 
and d = 0.49), (a) a randomly chosen initial condition; (b)-(g) 
observed steady states using the same initial condition in (a); 
(h) a distribution of steady states using four randomly chosen 
initial conditions (using 1000 runs each). 



(x, y, (/>)-space. As in ID cases, transient states will look 
like curved 2D surfaces mostly, but eventually the system 
will reach one of steady states, however long it takes [see 
Fig. [2|»]. When both the lattice and the opinion are 
periodic, non-flat metastable states can emerge, while 
there will be only flat steady states otherwise. When 
g < 1, opinions at a location can have more than one 
value: in other words, at a given location there will be 
a distribution of opinions. For non-flat steady states, 
the width of the opinion distributions at a given location 
will be finite, and as g gets smaller, the width of this 
distribution will increase. 



III. NUMERICAL RESULTS FOR THE 
PERIODIC CASE 

The toroidal (periodic) lattice with periodic opinions 
has richer dynamic behavior than other cases as we've 
seen in the previous section. In addition to the flat steady 
states, this system can have non-flat metastable states 
due to the periodicity of both the location and the opin- 
ion. We can categorize these non-flat states with period 
numbers with "period n^/n/" (n^,n; = 1,2,3,...; see 
Appendix for detail). For 2D cases, only one direction 
(either x or y) can be periodic. Even though any period- 
(n^/ni) steady states can exist, we only observed steady 
states of period [ni = 1) and period 1/ni (n^ = 1) 
mostly in our numerical simulations, and we will show 
them in later figures. For all numerical results except the 
one in Fig. [7^ a) , we assume that the density of agents is 
always one per location. 

In Fig.|3j we observe steady states that can emerge for 
TV = 1000, X = 1000, Y = 1, g = 1, and d = 0.49, with 
one randomly chosen initial condition (d = 0.49 is used 
as a maximal value instead of 0.5 because 0.5 is excluded 
as was discussed in the previous section). Since the dy- 
namics is stochastic, the system can reach the different 
types of steady states with the same initial condition. 
In addition to the flat steady states [Fig. [3j^b)] , non-flat 



FIG. 4. Scatter plots in two cases with 2D periodic lattices 
and periodic opinions, (a) TV = 10 000, X = 100, Y = 100, 
g = 1, and d = 0.49. (b) TV = 1000, X = 100, Y = 10, g = 1, 
and d = 0.49. 



steady states also appear [Fig. [3]jc)-(g)]. By repeating 
simulations using different sets of random numbers, we 
can find the distribution of types of steady states when 
we start from one given initial condition. The distribu- 
tions from different initial conditions don't have to be 
the same, as we've shown in Fig. [3](li) for four different 
randomly chosen initial conditions. 

How do we find out dynamical properties of the system 
with given parameters? Each initial condition will have 
its own distribution as we've seen in Fig. [3j therefore the 
distribution we obtain after averaging over those from 
all initial conditions characterizes the system. The more 
samples we choose and the more runs we perform for 
each sample, the more accurate this distribution should 
be. For each parameter set used subsequently, we will 
sample 1000 initial conditions randomly from the space 
of all initial conditions (in this case, {(xi, yi, < i < 
TV}), and run once each, to find out the approximate 
distribution (the total of 1000 runs). 

For 2D cases, the results are similar to those from ID 
cases. In Fig. |1J we looked at steady states for two differ- 
ent cases: (X, Y) = (100, 100) and (100, 10). In the first 
case, we observed two types of steady states: flat (~99%) 
and period 1 (~1%). While, in the second case, we ob- 
served four types of steady states: flat (~38%), period 1 
(-48%), period 2 (~12%), and period 3 (-2%). As we 
will show later in Fig. ^c), the width of the lattice in 
y-direction, Y , can change the dynamical behavior of the 
system, when X is fixed. When X = Y, the complete 
consensus was reached in almost all cases. 

Let us next observe how the threshold d changes the 
dynamic behavior of the system. In Fig. [5ja) , we com- 
pared our results with those from the Deffuant model 
by observing averaged numbers of opinion groups while 
varying d. Two agents are randomly picked at each 
time step, and g is set to 1. In numerical results with 
TV = 1000 (□), averaged over 1000 runs, gradual transi- 
tions were observed as the number of groups increases. 
We simulated for periodic opinions (o), too, and got al- 
most the same results. Two cases from our model were 
also simulated: a rectangular lattice with non-periodic 
opinions (O), and a toroidal lattice with periodic opin- 
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FIG. 5. (a) Number of opinion groups formed when varying 
d for four different settings: N = 1000, X = 100, Y = 10, 
g = 1, and d = 0.49, 1000 runs each. When counting numbers 
of groups, we counted groups with greater than 10% of the 
whole population, and ignored non-flat steady states, if they 
exist. The dashed line represents the prediction that, when 
d is in the range of l/(n + 1) < d < 1/n, there can be n 
opinion groups, (b) Distributions of observed steady states 
for the case with the toroidal lattice and periodic opinions, 
while varying d (N = 1000, X = 100, Y = 10, g = 1, and 
d = 0.49, 1000 runs each). 



ions (A) when N = 1000, X = 100, Y = 10, g = 1, and 
d = 0.49 for 1000 runs each. The results were similar to 
those from the Deffuant model, but numbers of groups 
tend to be a little smaller. 

In Fig. [5jb) , steady states for the case of the toroidal 
lattice with periodic opinions with N — 1000, X = 100, 
Y = 10, and g = 1, were observed when d is varied. There 
can be non-flat steady states. For flat steady states, we 
can categorize them with number of opinion groups: 1 
group, 2 groups, 3 groups, and so on. When d is less than 
~0.17, steady states are mostly flat, and as d gets smaller, 
the more opinion groups can exist as seen in Fig. [5ja). 
When d is between ~0.17 and ~0.27, fractional periodic 
states, mostly 1/2, emerge, while dominant steady states 
are those with 2 groups. When d is between ~0.27 and 
~0.32, this is where the transition occurs: the number of 
2-group steady states decreases quickly, while the number 
of period-1 steady states increases. When d is greater 
than ~0.32, all fiat steady states belong to the 1-group 
type, while there can be many types of non-flat steady 
states. Distributions don't change much as d increases 
up to 0.5. In the current case, steady states with period 
3 or higher don't appear much; however when Y/X is 
smaller, more types of periodic steady states will appear. 



2 groups 1 group period 1/2 
( c) d=0.35 (100 x 10, g= l) 



- ■ 



period 1 period 1 



1 group period 1 period 2 period 3 

FIG. 6. Scatter plots of some observed steady states for 
different d values from Fig.^b) (N = 1000, X = 100, Y = 10, 
g = 1, and d = 0.49, toroidal lattice with periodic opinions), 
(a) d = 0.15 (3 groups: 34%, 2 groups: 56%, 1 group: 2%, 
period 1/3: 6%, period 1/2: 2%), (b) d = 0.25 (2 groups 
52%, 1 group: 26%, period 1/2: 17%, period 1: 5%, period 1 
< 1%), (c) d = 0.35 (1 group: 41%, period 1: 48%, period 2 
10%, period 3: 1%). 



In Fig.|6j we show steady states in Fig.^b) for three d 
values: 0.15, 0.25, and 0.35.. When d is 0.15, three types 
of flat steady states dominate (~92%), while there exist 
non-flat steady states with fractional periods (mostly 1/2 
and 1/3). When d is 0.25, more than 3/4 of steady states 
are flat (~78%) still, while non-flat steady states with 
period 1/2 and 1 also exist. Here there are two kinds 
of period-1 steady states: one has one band (n^ = 1 
and ni = 1), and the other has two bands (n^ = 2 and 
ri; = 2), which is possible because d is small. When there 
are two bands, there are two disjoint groups of agents, 
even though each group has agents with a full spectrum of 
opinions, and they were observed only in the approximate 
range of d between 0.18 and 0.25. When d is 0.35, only 
1-group flat states were observed, and non-flat steady 
states of periods 1, 2, and 3 were also observed. 

Finally, we can ask how other parameters will influ- 
ence the outcome. In Fig. [TJa), N varies from 100 to 
2000 when X = 100, Y = 10, g = 1, and d = 0.49 
to observe how the density [N/ (XY)] changes dynamical 
behavior. When the density is much smaller than 1, the 
number of non-flat steady states decreases. But we can 
clearly observe that if the density is greater than about 
0.5, distributions do not change much. This result shows 
that the density doesn't have to be high to find out dis- 
tributions of steady states, and that the density of 1 is 
good enough. 

In Fig. ^b), we vary the size of the lattice while keep- 
ing the density and the shape fixed. As X increases from 
100 to 500 when N = XY, Y = X/10, g = 1, and 
d = 0.49, the distributions do not change much. We can 
interpret this result as a sign that there is no size effect 
unless the size is too small. 

In Fig. [7](c), Y varies from 10 to 300, changing the 
shape of the lattice, when N = XY, X = 100, g = 1, 
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FIG. 7. Distributions of steady states varying TV, X, Y, 
and g (d is fixed at 0.49, using 1000 runs for each case), (a) 
Varying N (X = 100, Y = 10, g = 1). The density changes, 
(b) Varying TV, X and V (TV = XY, Y = X/10, g = 1). The 
density and YjX do not change, and the size effect can be 
seen, (c) Varying TV and Y (TV = XY, X = 100, 5 = 1). The 
density and X do not change, while the shape of the lattice 
changes, (d) Varying g (TV = 1000, X = 1000, Y = 1). 



and d = 0.49. As Y approaches the value of X, the num- 
ber of non-flat steady states decreases. As we've seen in 
Figs. [4] and [6j the periodic behavior only appears in one 
direction in the case of toroidal shapes. If opinions are 
periodic in x, opinions along the y direction for a given x 
value are more or less constant. When Y is 10, periodic 
behavior can be only seen in x, but as Y increases, non- 
flat steady states along the y direction starts to emerge. 
When X — Y, steady states can be periodic either in x 
or in y with equal probabilities. When Y is greater than 
X, the likelihood of forming steady states that are peri- 
odic in y will be greater, and the probability of getting 
the periodic steady states starts to increase again. The 
more elongated the shape of the lattice, the greater is the 
likelihood of finding periodic steady states. 

In Fig.[7^d), the convergence parameter g varies when 
TV = 1000, X = 1000, Y = 1, and d = 0.49. As g gets 
smaller, the number of non-flat steady states decreases. 
To make those periodic steady states disappear, g has 
to be very small, ~10~ 4 , in this case. The convergence 
parameter controls how fast opinions converge; in ad- 
dition, when g is small, the distribution of opinions at 
the same location for non-flat steady states gets wider 
because agents can move farther away from a location 
without changing their opinions much. 

In Fig. U we observe widths of distributions of opin- 
ions at a given x for some of non-flat steady states from 
Fig.fjjd). In Figs. ||a)-(g), non-flat steady states when 
q = 10 _1 . 10 -2 , 10 , and 10~ 4 are shown. Insets show 



distributions of opinions <f> at x = 500. We use opinions 
of agents residing at x = 500 at a number of different 
time steps after steady states are reached (for example, 
1000 time steps with 100 time steps apart) [55] • If these 
frequencies are fitted to Gaussian functions as shown in 
insets in Fig.|8j the standard deviations, a, will represent 
the widths of these distributions. We can argue that a 
is a measure of stability. As discussed in Sec. II, we can 
call non-flat steady states metastable, because lifetimes 
of structures formed in (a;, y, 0)-space are finite, but ex- 
tremely long. In other words, the probability of decay 
within a certain finite time period is extremely small, 
even though it will not be zero exactly as in flat steady 
states. If we define R as the ratio of the width of a band 
to the distance between two adjacent bands along the x 
direction, 



R - — 777 — n ' CT ' 



(7) 



where n/A/n^ is the inverse of the slope of the line in 
(x,4>). Since ni = 1 always in current examples, a = R. 
If R <?C 0.5, the structure is quite stable. If R > ~0.5, ad- 
jacent bands will overlap and will be broken immediately. 
If R comes close to 0.5, the state is likely to become flat or 
morph into another type of steady states in a short time 
period. In general, the bigger a is, the more unstable the 
state is and the less likely it will be reached. 

How is a affected by parameters? We have observed 
111 Fig. [8j that the greater g, the smaller is a. This width 
also depends on the period number n^/nf, for example, 
the width of opinions of period-2 steady states is twice as 
wide as that of period-1 steady states, because the opin- 
ion difference between two adjacent locations is twice as 
big for period-2 steady states. Then we can generalize 
that for the same g value, the width of opinions at a lo- 
cation for steady states with period / rt\ is proportional 
to 7i0 / ni . This explains why steady states with higher pe- 
riods disappear quickly as g gets smaller. When g = 1, a 
comes close to zero unless N is very small. That leads us 
to claim that N can play a role, too: as N gets smaller, 
a becomes wider because, when N is small compared to 
XY, the possibility of having empty locations increases 
and an agent can move further away without changing 
its opinion [see Fig. |7|[a)]. In addition, stability of non- 
flat states also depends on YjX [see Fig. |7jc)]. We can 
summarize based on our results so far that, in general, 
the stability of a certain periodic state depends directly 
on Nj '(XY), YjX, g, and n/, while d determines what 
types of periodic states are allowed to exist. 

Figure [8jh) shows a special case. It has three opin- 
ion groups, but unlike opinion groups found in flat states 
when d is small, they have non-zero widths of opinions 
and an opinion of each agent is not stationary. In addi- 
tion, the lifetimes of these states are not long. In most 
cases, this type of states stayed intact for more than 10 
times of that it took for most steady states to be formed. 
This type of states were observed only in the case with 
g = 10~ 4 for about one percent [note that since they 
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(a) g^lO'^er.l, t=5xlO s ) (b) g=10' 2 (per.l, t=5x!0 s ) (c) g=10' 2 (per.2, t=5x!0 5 ) (d) g=10~1(per.3,t=5xi0 5 ) 






c>=0.022 




(e) g=10 (per.l, t=10 6 ) (f) g=10' 3 (per.2, t=10 6 ) (g) g^lO^Cper.l, t=10 6 ) (h) g=10 (3 groups, t=10 ) 






ct=0.012 
^Af ■*'''<***'' -0.016 

A A A : 



500 



500 * 



500 x 



FIG. 8. Scatter plots of non-flat steady states from Fig. [7|d). The widths were obtained by fitting the distributions of 4> a t 
x = 500 (pointed by arrows) by Gaussian functions, and a is the standard deviation (shown in insets). 



eventually converge to flat states, they were counted as 
flat states in Fig. j7^d)] . They could emerge, because pe- 
riodic opinions are used and opinions from three groups 
can be balanced in a circle. Also, g has to be small 
enough; if g is bigger, widths of opinions will be bigger 
and the structure becomes unstable. 



IV. DISCUSSIONS 

We introduced a simple model of opinion dynamics 
that uses only nearest-neighbor interactions. To repre- 
sent the geometrical space we live in, we used 2D lattices, 
on which agents move randomly. Both the lattice and the 
opinion can be periodic or non-periodic. We explored 
some regions of the parameter space, and found rich dy- 
namic behavior, especially when both are periodic. One 
might argue that periodic opinions and toroidal lattices 
are not realistic and even artificial. But the opposite 
might be true: linear opinions and rectangular lattices 
are rather special cases. The surface of the Earth is finite 
and has no boundary. Opinions, discrete or continuous, 
are not always linear, either. In some cases, opinions can 
be better represented by more general structures other 
than a line. In short, our results show that if we gener- 
alize spaces for opinions and lattices, dynamic behavior 
of those systems can be richer. 

In reality, a consensus is not reached easily. Con- 
sider political opinions as an example. Even though 
political parties have been formed in advanced soci- 
eties, people in one party usually have various politi- 
cal views. Another example is the existence of dialects, 
if we regard languages as opinions. These phenomena 



can be explained by inherent heterogeneity of agents, 
or bounded confidence. Our model, however, adds an- 
other explanation: metastability through local interac- 
tions. Metastable states in our model indicates that 
locally-converged, but not globally-converged, states can 
be sustained for a long time in certain situations. For 
example, if we are surrounded by like-minded people, we 
seldom change our views and even believe that everybody 
is similar to us, which isn't true in general. 

In general, steady states or equilibrium states in mod- 
els that are closed are hard to be realized in real sys- 
tems like societies because societies are fundamentally 
open and extremely noisy. However, behavior of tran- 
sient states and emergence of different types of steady 
states can shed some light on understanding how real 
systems behave. In this model, metastable states, which 
have the full spectrum of opinions, are only observed in 
cases of toroidal lattices with periodic opinions. But even 
for cases with rectangular lattices, as the size of the lat- 
tice gets bigger, the overall behavior of locally-converged 
transient states seen in Fig. [2] will be similar as long as 
the opinions are periodic. Then we can interpret our 
real-life states containing the wide spectrum of opinions 
as transient states that are moving slowly toward steady 
states. 

The model considered here can be expanded or modi- 
fied. We can assume g and d are not constant throughout 
the whole population (so called, heterogeneous agents [T2l 
Hl?l |2"5] ) . The definition of neighbors can be modified by 
using a bigger range, so that the network structure be- 
comes more realistic. We can also consider an additional 
co-evolving network like an acquaintance network on top 
of our model, and it will be investigated in a future arti- 
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1/3 1/3 1/3 

1/3 A 1/3 1/3 A 1/3 1/3 A 1/3 

^~iir* *~iir* *~ilr* 

position: x—\ x x+\ 

opinion: <f(jc — 1) §(x) §(x+l) 

FIG. 9. Schematic diagram for an ID case (Y — 1) with 
g = 1 and the maximal d. 

cle. 

The models for social dynamics like this are not devised 
for predicting the future of specific systems, but can help 
us understand dynamical properties. In addition, they 
can be adapted to a wide variety of different systems 
with similar features. 



ACKNOWLEDGMENTS 

The author thanks Chang- Yong Lee for his help on 
finding out that some unexpected results were actually 
metastable states, while at first it had seemed to be from 
unknown numerical errors. The author also thanks L. 
E. Reichl for her thoughtful comments concerning this 
manuscript. This work was supported by Kongju Na- 
tional University. 



Appendix: Solving an ID case 

We make the lattice one-dimensional by setting Y = 
1, and we also assume g = 1 and the maximal d (d — 
1 for non-periodic opinions, and d = 0.49 for periodic 



opinions). If the density of agents, the average number of 
agents at each location, is high enough, we can simplify 
the system by making the opinion to be the function 
of the location, ignoring the states of individual agents. 
That is due to the fact that every agent in the same 
location will have the same opinion instantly. Then the 
overall Y-agent dynamics is reduced to the dynamics of 
</>(#), an Y-dimensional vector. 

Since an agent has three choices to move (left, right, 
and staying) with the probability 1/3 for each move (see 
Fig. [9]), the opinion at x, (f){x), at the next time step will 
become [4>(x — 1) + (f>(x) + (f>(x + l)]/3. This is a discrete- 
time linear dynamical system, which can be character- 
ized by the tridiagonal matrix with all non-zero elements 
1/3. Usually we are interested in finding states where 
4>(x) does not change with time, and they are eigen- 
states of this matrix with the eigenvalue 1. If we look 
at the dynamics locally, when (f>(x — 1) = 4>{x) — 5 and 
(f)(x + 1) = (j)(x) + S, (j)(x) will not change at the next 
time step, where 6 is a small real number at x. 

Since both x and </> are bounded, we cannot find eigen- 
states easily, except one trivial case when <p(x) is constant 
with respect to x (5 = for all x). If a system arrives at 
these states, a consensus has been reached. In most cases, 
these are the only eigenstates, but when both the lattice 
and the opinion are periodic, the periodic eigenstates can 
exist as long as the boundary conditions are met, where 
S is a non-zero constant for all x. However, S needs to be 
certain discrete values, and we can name these periodic 
eigenstates with period numbers. The "period n^/n" 
(n^, ni — 1,2,3,.. .) means that opinions have cycles 
of changes while there are n; cycles of changes along the 
x direction; hence only when S = {n^/ni)X, they will 
become eigenstates. 
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